Astronomy & Astrophysics manuscript no. PDM' FINAL 


©ESO 2013 


January 14, 2013 





Signatures of radial migration in barred galaxies: 
Azimuthal variations in the metallicity distribution of old stars 

P. Di Matteo 1 , M. Haywood 1 ^. Combes 2 , B. Semelin 2 , O. N. Snaith 1 



CO 

T— I 

o 



< 

o 

6 

c3 



> 
in 

T— I 

o 

m 



X 



1 GEPI, Observatoire de Paris, CNRS, Universite Paris Diderot, 5 place Jules Janssen, 92190 Meudon, France 
e-mail: paola . dimatteo@obspm . f r 

2 LERMA, Observatoire de Paris, CNRS, 61 Av. de 1' Observatoire, 75014 Paris, France; Universite Pierre et Marie Curie, 4 place 
Jussieu, 75005 Paris, France 



Accepted, Received 



ABSTRACT 



By means of N-body simulations, we show that radial migration in galaxy disks, induced by bar and spiral arms, leads to significant 
azimuthal variations in the metallicity distribution of old stars at a given distance from the galaxy center Metals do not show an 
axisymmetric distribution during phases of strong migration. Azimuthal variations are visible during the whole phase of strong bar 
phase, and tend to disappear as the effect of radial migration diminishes, together with a reduction in the bar strength. These results 
suggest that the presence of inhomogeneities in the metallicity distribution of old stars in a galaxy disk can be a probe of ongoing 
strong migration. Such signatures may be detected in the Milky Way by Gaia (and complementary spectroscopic data), as well as in 
external galaxies, by IFU surveys like CALIFA and ATLAS3D. Mixing - defined as the tendency toward a homogeneous, azimuthally 
symmetric, stellar distribution in the disk - and migration turns out to be two distinct processes, the effects of mixing starting to be 
visible when strong migration is over. 

Key words. Galaxies: abundances; Galaxies: evolution; Galaxies: structure; Galaxies: kinematics and dynamics; Methods: numerical 



1. Introduction 

Stellar asymmetries in galaxy disks, such as bars and spiral pat- 
terns, have been known for a long time to drive galaxy evo- 
lution. Stellar bars, in particular, constitute an efficient way 
to redistribute angular momentum amon g the different galaxy 
components - gas, stars and dark matter (iBournaud & Combesl 
l2002t iBerentzen etaD l2007h lAthanassoula & Misiriotisl l2002h . 
Angular momentum redistribution is particularly efficient at res- 
onances, and a coupling between the action of a bar and that 
of sp iral arms can cause a rapid migration of stars through the 
disk dMinchev & Famaevl201Ql:lMinchev et al.ll201 ll) . with char- 
acteristic time scales significantly lower than those predicted 
when the action of spiral arms alone is taken into account 
(ISellwood & Binnevl |2003: iRoskar et all l2008al) . Mergers and 
perturbations from sa tellite galaxie s can a l so induce radial mi- 
gration , as shown by lOuillen et al.l (I2009I) . iBekki & Tsujimotol 
(1201 ll) and lBirdetal.1 (120121) . 

If radial migration has received considerable attention in 
the last few years, driven by the need to interpr et a num- 
ber of galactic and extragalactic observations ([H aywood 
20081: iBovv et al.ll2012l: iB^kos et al.ll201 l l : lYoachim et al.ll2012l: 
Radburn-Smith et al.ll2012l) . it is still not clear how important 
this process is in galaxy evolution, whether all disk galax- 
ies go through a phase of strong mi gration, and how th e 
Hubble type and bar strength affect it. iBrunetti et al.l (1201 ll) . 
for example, showed that not all barred galaxies experi- 
ence strong diffusion, this depends on the bar strength, and 
thus ultimately on the stability of the disk. If more work 
is needed to quantify systematically the importance migra- 
tion has in galaxy evolution, it is also essential to understand 
the signatures this process leaves in a galaxy. Truncated stel- 
lar density profiles, upturns in the stellar age profiles, flat- 



tening of stellar metallicity gradients, thick disks, et c. 
be the result of migration processes (iDebattista et al 
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Roskar et al. 2008b: lLoebman et al.ll201 ll iMinchev et al.l 



2006 



2011 



20121) . Unf ortunately, however, they are not uniquely caused by 
migration (lElmegreen & Hunterl 1200 6: Sanc hez-Blazquez et al.l 
120091: lOu et al.ll201 ll IBournaud et al.ll2009l) . Is there any other 
way to quantify the strength of the process, and to understand if 
a galaxy is currently experiencing a strong redistribution of stars 
in its disk? 

In this work we report on the results of a numerical study 
which shows that in the phase of strong migration induced by 
a bar galaxy disks exhibit significant azimuthal variations in 
the metal distribution of their old stellar component. These az- 
imuthal variations remain strong during the phase of high bar 
strength (~ lGyr in our simulations), when radial migration is 
maximum, and decrease later on. At this point, the inner stellar 
disk becomes homogeneous, on kpc scales, within few rotational 
periods after the end of the migration phase. We suggest that the 
presence of azimuthal inhomogeneities in the metal distribution 
of old stars in a barred galaxy can be used as evidence that its 
disk is going through a phase of significant stellar migration, re- 
lated to a strong bar. 

The paper is organized as follows: in Sect]2]a description of the 
initial conditions and numerical method adopted for the simu- 
lation are given; in Sectj3]the results are presented, firstly dis- 
cussing the mode and timing of radial migration in barred galax- 
ies, and then presenting the effect of this migration on the redis- 
tribution of metals in a stellar disk, together with a discussion on 
its dependency on the initial metallicity gradient; in Sect|4] the 
main conclusions of this work are described. 
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2. Initial conditions and numerical method 

We report on the results of a collisionless N-body simulation of 
an isolated galaxy, with an initial bulge-to-disk ratio equal to 0. 1 , 
and no gas in its disk. This simulation is one of a set of three dis- 
sipationless high-resolution simulations, with varying bulge-to- 
disk ratio. We present the analysis only for one of the three sim- 
ulations, since the results and conclusions found are common to 
all the cases analyzed. Note that hydrodynamics is not included 
in this work, since we are interested in quantifying the effect of 
radial migration alone. Studying the complex interplay between 
gas, induced star formation, and stellar dynamics is beyond the 
scope of this paper. However, for a Milky Way-like galaxy, with 
a quiescent star formation of few M©/yr, we do not expect the in- 
clusion of star formation to drastically change the results, since 
the amount of newly formed stars per rotation constitutes only a 
small percentage of the total stellar mass in the disk. 

The halo and the bulge of the galaxy are modeled as a 
Plummer sphere (iBinney & Tremainelll987l) . with characteris- 
tic masses M B = 9 x 10 9 M o and M H = 1.02 x K^AfcQ and 
characteristic radii re = 1.3 kpc and rH=10 kpc. The stellar disk 
follow s a Miyamoto-Nagai density profile (iBinney & Tremainel 
fT987h . with mass M* = 9 x 1O 1O M and vertical and radial scale 
lengths given, respectively, by h*=0.5 kpc and a* =4 kpc. The 
initial disk size is 13 kpc and the Toomre parameter is set equal 
to Q=1.8. The galaxy is represented by N tot = 30720000 parti- 
cles redistributed among dark matter (Nu = 10240000) and stars 
(Altars = 20480000). To initialize particle velocities, we adopted 
the method described in lHernquistl (Il993l) . 

To model the galaxy evolution, w e empl oyed a Tree-SPH 
code, presented in ISemelin & Combesl (120021) and we refer the 
reader to this paper for a full description. For the dissipation- 
less simulations analyzed in this paper, the SPH part of the code 
has been switched off, and gravitational forces are calculated us- 
ing a tolerance parameter 6 = 0.7, including terms up to the 
quadrupole order in the multiple expansion. A Plummer poten- 
tial is used to soften gravitational forces, with a constant smooth- 
ing length 6 = 50 pc for each of the different species of particles. 
With this spatial resolution, it is possible to resolve the vertical 
structure of thin disks, and follow small scale inhomogeneities in 
the disk. The equations of motion are integrated using a leapfrog 
algorithm with a fixed time step of At = 0.25 Myr. 

Stars in the disk are assigned a metallicity z m which de- 
pends on their initial distance R from the galaxy center, accord- 
ing to the formula z m = zolO~°' 07R , zo - 3z© being the central 
(R = 0) metallicity, and Zq the solar value. We thus assume ini- 
tially no azimuthal variation at any given radius. For more details 
on this choice, see Appendix lAlThe initial metallicity gradient 
has been also varied from its value A\ Fe i m = -0.07 dex/kpc 
dMaciel & Koppenlll994l : Ivan Zee etal]ll998h to study the de- 
pendence of the results on this initial choice, as discussed in 
SectES 

3. Results 

3.1. Radial migration in barred galaxies: how and when? 

Due to the high number of particles employed and the re- 
duced noise in the gravitational force field, it takes about 
0.8 Gyr to see the appearance of asymmetric structures in 
the stellar disk. Indeed, only at this epoch strong m = 2, 4 

1 This value corresponds to the mass inside a sp here of radius R= 
35 kp c, and guarantees realistic rotation curves (see IChilingarian et all 
1201 Ol for a similar set of galaxy models) 
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Fig. 1. Strength of the A\ (black line), A 2 (red line) and A 4 (blue 
line) stellar asymmetries as a function of time. 
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Fig. 2. Density maps of the stellar disk, seen face-on, at different 
times. Each panel is 70 kpcx70 kpc in size. 



asymmetries, associated with a bar and spiral structure, develop, 
maintaining a nearly constant strength for about 1 Gyr (see 
FigQ] and [2]). The decline that follows in the bar and spiral 
arms strength at about t=2 Gyr is associated to a conspicuous 
vertical buckling of the bar, and the subsequent formation of a 
boxy/peanut shaped bulge, as obs erved many times in N-body 
simulations dCombes & Sandersl [l98lt ICombes et alJ 11990: 
iMartinez-Valpuesta et al.ll2006UAthanassoulall2008l) . It is worth 
noting that at later times the galaxy is still barred, but with 
isophotes rounder than those characterizing the phase of high 
bar strength (Fig. 0). 

Motivated by the results of Fig. [T]|2l in the following, when 
discussing the disk evolution, we will refer to some distinct 
temporal phases: early phase (t < 0.8 Gyr); bar growth (very 
rapid phase, from t = 0.8 Gyr to t = 1 Gyr); strong bar activity 
(corresponding to the time of high, nearly constant bar strength, 
i.e. 1 Gyr< t < 2 Gyr); weak bar (t>2 Gyr). 
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Fig. 3. Face-on density distribution of stars initially in a radial 
annulus 7 kpc< R <9 kpc ( in each panel, the average value, 
R = 8 kpc, is indicated, by a dashed white circle). 



At the epoch of bar formation, and during the whole 
phase of strong bar activity, stars in the disk experience a 
significant spatial redistribution, with a probability of migration 
m aximal at the bar res o nance s, in agreement with the findings 
of iMinchev & Famaeyl (I2010I) . Fig. [3l for example, shows the 
spatial distribution over time of stars that at the beginning of 
the simulation were in the radial annulus 7 kpc < R < 9 kpc, R 
being the distance from the galaxy center. This region includes 
the corotation radius at early times (the bar corotation is 
initially located at R ~ 8 kpc, but it moves rapidly outwards 
to R ~ 10 kpc). Before asymmetric structures appear in the 
disk, stars initially belonging to this annulus maintain a radially 
symmetric distribution in the plane. As soon as the bar starts to 
develop, the spatial distribution of stars becomes asymmetric: 
elongated along the bar axis in the inner regions, and showing 
signs of spiral structures in the outer regions. Stars initially 
confined in this annulus can migrate toward the inner and outer 
disk. In particular, stars moving outwards do so via spiral arms, 
which form together with the bar. These arms are visible during 
the phase of strong bar activity, fading away only 2 Gyr after 
the beginning of the simulation, when the bar strength declines 
(note that the length of the bar is always < 9 kpc, as it can be 
deduced by the location of the two density maxima -yellow 
clumps in Fig. [3]for t > 1.1 Gyr - which identify its end). 
As we will discuss in the next sections, the way radial migration 
occurs, through these asymmetric structures, has an important 
impact on the spatial redistribution of metals and its variation 
with time. 

Fig. [3] clearly shows that in these simulations radial migra- 
tion is initiated by the growth of the bar, at t ~0.8 Gyr. But 
for how long does it last? Is radial migration only a transient 
phenomenon, related to the formation of the bar, or does it 
last longer than few hundreds of Myrs? To answer this ques- 
tion, we have evaluated the fraction of stars that migrate of 
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Fig. 4. Fraction of stars that migrate of | AR \> 3,5,7 kpc from 
their location R at the beginning of the time window, at four 
different time intervals (from top to bottom): t e [0. - 0.4] Gyr; 
t g [0.4 - 1.4] Gyr; t e [1.4 - 2.4] Gyr; t e [2.9 - 3.9] Gyr. Bins 
are 0.1 kpc in width. 



more than | AR \=\ R(t f ) - R(ti) \= 3,5,7 kpc from their 
location at the beginning of the time window [ti,tf], with 
[t i9 tf] = [0., 0.4], [0.4, 1.4], [1.4, 2.4] and [2.9,3.9] Gyr, respec- 
tively (Fig. [4]). Before asymmetries appear in the stellar disk, the 
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fraction of strong migrators (| AR \ > 5 kpc) is null all over the 
disk, except at the outer edge, where an initial relaxation from 
the initial conditions leads the disk to slightly expand with re- 
spect to the truncation radius initially at R = 13 kpc. The frac- 
tion of migrators drastically changes in the following Gyr, from 
t = 0.4 to t = 1.4 Gyr: the whole outer disk is affected by a 
significant redistribution of stars, with up to 30% of stars in this 
phase moving by more than 7 kpc. There are two radii where 
the relative fraction of migrators is higher: at R = 8.5 kpc and 
R = 12.5 kpc, corresponding respectively, to the corotation and 
outer Lindblad resonance at that time. If migration is strong at 
the epoch of bar growth, it is still significant at later times: in the 
time interval t - 1 .4-2.4 Gyr, mostly corresponding to the phase 
of strong activity of the bar, 20% of stars migrate by more than 
| AR \> 1 kpc from the corotation (located at R = 10 kpc at this 
epoch). A weaker, but still not negligible migration, is present 
also at later times (t = 2.9 - 3.9 Gyr), when the bar strength has 
diminished, after the vertical buckling. 

We note that the distribution of the fraction of migra tors versus 
radius is significantly different from that predicted by lBird et al.l 
(120121) (see Fig. 4 in that paper). While in their isolated galaxy 
models, this fraction follows the surface density profile of the 
disk, decreasing with radius, our models predict the opposite 
trend, with a very low probability of strong migration in the in- 
ner disk - where the bar traps a large fraction of stars - and local 
maxima in the outer disk associated to the location of the bar 
resonances. This discrepancy is probably related to the fact that 
their isolated disks are stable against bar formation, and, as a 
consequence, the effect of scattering at the bar resonances is not 
taken into account. 



3.1 .1 . Comparison with Minchev et al results 

The results discussed in Sect. I3.ll agree with those already pre- 
sented in a number of previous works, concerning the study 
of radial migration in b arred galaxies (iBrunetti et"aD l201ll: 
iMinchev & Famaey 20]1; IMinchev et all 1201 ll l2012jT In par- 
ticular, as discussed in several works by Minchev and collabora- 
tors, we can confirm that radial migration in these simulations is 
not only associated to the epoch of bar growth: migration is in- 
deed initiated at that time, but it lasts for all the phase of strong 
bar activity, and also, at lower levels, later on (for t > 2Gyr). 
Moreover, the stability of the simulated disks over several hun- 
dreds of Myrs (see Fig. [T] and 13 guarantees that radial migra- 
tion observed in these simulations is not a consequence of ini- 
tial violent instabil ities in the disk, as recently suggested by 
Roskar et al. d2012h to be the case for the simulations analyzed 
bv lMinchev et al.1 (1201 l[l2012h . 



3.2. Azimuthal variations in the metallicity distribution 

In the previous subsection we have seen that the presence of 
asymmetries in the stellar disks causes a significant radial redis- 
tribution of stars. This redistribution lasts all the time the bar is 
strong, and declines to lower levels after the bar buckling. Since 
migration is efficient mostly in the phase of strong bar activ- 
ity (Fig. [4]), is there any observational signature that can help in 
quantifying the strength of the bar? In other words, is there a 
way to quantify if the galaxy is experiencing a strong phase of 
bar activity and thus of strong migration? 
In this section and in the following, we will show that the pres- 
ence of azimuthal variations in the metallicity distribution of old 




5 10 15 

R [kpc] 



Fig. 5. Azimuthally averaged metallicity profile of the stellar 
disk at different times, as shown in the legend. The initial metal- 
licity gradient is A[p e /H] = -0.07 dex/kpc. Metallicity disper- 
sions at three different times (see corresponding colors), and at 
two different radii (R = 0, thin bars; R = 8 kpc, thick bars) are 
shown in the top-right part of the plot. 
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Fig. 6. 2D [Fe/H] maps of the stellar disk, seen face-on, at dif- 
ferent times. Each panel is 70 kpcx70 kpc in size. 



stars in a disk is a powerful probe of the state of bar activity and 
of its impact on the stellar disk. 

We will start this analysis by discussing the evolution with 
time of the stellar metallicity profile (see Fig. [5]). In the phase 
of strong migration, the initial metallicity gradient flattens out 
significantly in the outer disk regions, varying from its initial 
value A[Fe /Hi = -0.07 dex/k p c to ArFe/Hi = -0.01 dex/kpc 
(see also iFriedli et al] (11994 : IMinchev etaD (12011b ). while 
changing only marginally in the inner regions (from A[ Fe /H] = 
-0.07 dex/kpc to A[p e /H] = -0.06 dex/kpc). This leads to a char- 
acteristic change of slope in the metallicity p rofile, as often ob- 
serve d in disk galaxies (for a recent work, see lScarano & Lepinel 
120121) . At the same time, the metallicity dispersion increases. At 
t=0, by construction, the dispersion is null; at t=0.4 Gyr, well 
before bar formation, the dispersion around the mean metallicity 
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Fig. 7. Maps of azimuthal metallicity variations, £[F e /H] (dex), 
at the different times shown also in Figj6] Black contours cor- 
respond to equally spaced isosurfaces of the differential stellar 
density, 



caused by epicycle "blurring" is about 0.12 dex0 and as soon as 
the bar forms, the value raises up to about 0.16 dex, and stays 
constant till the end of the simulation. Radial migration pro- 
cesses related to the phase of strong activity of the bar thus add 
an extra dispersion with respect to the pre-migration value, of 
about 0.1 dexB 

In Fig(5] we have evaluated the azimuthally averaged metallic- 
ity profile. But how far from axisymmetry is the real metallic- 
ity distribution? In other words, what is the variation around the 
mean value at any given distance from the galaxy center? The 
face-on maps of the stellar [Fe/H] distribution at different times 
are shown in Figj6] Once formed (between 0.8 and 1 Gyr) the 
bar is clearly recognizable in the inner regions, where the metal- 
licity, initially axisymmetric, becomes elongated along the bar 
major axis. Also, the distribution in the outer regions is far from 
axisymmetry, showing signs of spiral structure and metallicity 
inhomogeneities up to the edge of the disk. 

These inhomogeneities are related to the way radial migra- 
tion occurs in galaxies, through spiral patter ns (see Sect. 13.11 
and iMinchev et all 120121: [Grand et al.l l2012h : metal-rich stars 
which move to the outer disk are m ostly from the region outside 
the corotation (iBrunetti et al.l 1201 Ik and migrate through spiral 
patterns to the outer parts of the disk. In other words, migration 
is not axisymmetric. To show that the variations observed 
in the metallicity maps are related to the appearance of the 
bar and spiral arms, and that they fade away with the stellar 
asymmetries, we present maps of the azimuthal metallicity 
variations, £[p e /H], evaluated at different times during the galaxy 

2 The value is measured at R = 8 kpc, but note that it depends on the 
distance from the galaxy center. At R = 0, for example, the metallicity 
dispersion in equal to 0.06 dex at t=0.4 Gyr and changes only slightly 
at the epoch of strong bar activity, increasing up to 0.08 dex 

3 This value has been obtained by quadratic sum, since the two dis- 
persions, before and after migration, are independent. 



evolution (Fig|7]). Here tf[F e /H] = [Fe/H]- [Fe/H] axy , with [Fe/H] 
the raw value of the metallicity and [Fe/H] ax}; the azimuthal 
average, evaluated, at each time, in radial bin of 0.1 kpc. These 
variations are compared with the stellar density distribution 
at the same epoch, quantified through the differential stellar 
density £^yy=(2 - H axy )/I< aX y, where £ is the raw value of the 
stellar density and 2 ^ is the azimuthally averaged density at 
each time (see also IMinchev et al] 120121) . These plots show 
that azimuthal variations in the metallicity of the old stellar 
component appear as soon as stellar asymmetries start to de- 
velop. As the bar develops, £[F e /H] maps become quite complex, 
with maxima (and minima) inside the bar radius aligned with 
(mostly perpendicular to) the bar major axis. Also, the outer 
disk - which has expanded in size since the epoch of bar 
formation, in agreement w ith previous works (among others, 
IMinchev etaD 1201 lL l2012h - is characterized by metallicity 
variations, the absolute magnitude of which can reach 0.1 dex for 
an initial gradient A[p e /H] = -0.07 dex/kpc. Strong metallicity 
variations are present in the inner disk for all the time that 
stellar asymmetries are strong, and they fade away as the bar 
and spiral arms strength diminishes. However, while in the bar 
region metallicity variations are always aligned with the bar, in 
the outer disk maxima and minima do not necessarily coincide 
with maxima and minima in the differential stellar density. This 
behavior is due to the fact that once they have migrated into 
the outer disk stars have an angular rotation velocity different 
from the pattern speed of these asymmetric structures, and can 
thus intercept the arms at different location during their motion 
around the galaxy. 



3.3. Dependence on the initial metallicity gradient 

The amplitude of the inhomogeneities, formerly discussed, de- 
pends on the initial metallicity gradient of stars in the disk, 
as shown in Fig|8] (top panel), where the value of £[Fe/H] is 
shown as a function of azimuth, for three different values of 
the initial gradient, steep (A[p e /H] = -0.19 dex/kpc), moder- 
ately steep (A[Fe/H] = -0.07 dex/kpc) and nearly flat (A[p e /H] = 
-O.Oldex/kpc). At t=l.l Gyr, when the bar is strong, the maxi- 
mum variations (of the order of 0.2 dex) are found for the case 
of steep gradient, while very small variations (~ 0.01 dex) are 
produced if the initial gradient is nearly flat. Once the phase of 
strong activity of the bar is over (t=3.95 Gyr), the disk is chemi- 
cally homogeneous, with the metallicity variations fading toward 
zero for any choice of the initial gradient. 
It is interesting to note that a small £[Fe/H] does not necessarily 
imply a weak or absent radial migration. Small variations, due 
to an initial nearly flat metallicity gradient, may indeed hide a 
strong migration. 

To quantify the amount of spatial redistribution of stars 
in a disk it is indeed necessary to normalize to the value of 
the gradient, as shown in Fig [8] (bottom panel). Here the ratio 
£[F<?/#]/A[Fe/H] is plotted as a function of time, for three different 
values of the initial metallicity gradienfl In evaluating this ratio, 
we have used the initial value of the metallicity gradient, A[p e /H] » 
since it varies only marginally over time in the bar region^ (cf 
FigI5]and discussion in Sect. 13.21) . 



4 In this plot, 6[Fe/H] represents the maximum absolute value (over 6) 
of the abundance variation at any given radius R 

5 Note that this choice leads to underestimate the ratio S[ Fe / H] / A [¥e /n] 
in the outer disk, where the gradient flattens considerably after bar for- 
mation. 
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Fig. 8. Top panel: Metallicity variations, £[Fe/H], as a function of 
the azimuth, for three different values of the initial metallicity 
gradient: A[p e /H] = -0.01 dex (blue line), A[p e /H] = -0.07 dex 
(red line), A[p e /H] = -0.1 dex (green line), and at four differ- 
ent radii in the galaxy disk. Variations are shown at the time 
of strong bar activity. Note that the values corresponding to 
R = 7 kpc, R = 9 kpc, R = 12 kpc have been arbitrarily shifted of 
0.5 dex, 1 dex, 1.5 dex with respect to their real values. Bottom 
panel: £[F e /H]/A[F e /H] ratio as a function of time, for an initial 
metallicity gradient A[p e /H] = -0.01 dex/kpc. The ratio is shown 
at three different disk radii: R = 3 kpc (black line), R = 9 kpc 
(green line), and R = 12 kpc (blue line). Note that this ratio is 
independent of the gradient, as detailed in the text. 



As expected, Fig [8] shows that this ratio is independent on the 
initial metallicity gradient, and it is a tracer of the strength of 
radial migration, i.e. of the contamination, at any given radius, 
of stars coming from different disk radii. As a first approxima- 
tion, this value can be seen as the difference between R and R, 
with R being the radius where migrating stars come from, i.e.: 
£[Fe/H]/A [Fe /H] = (A [Fe /H]^ - A [Fe /H]^)/A[Fe/H] = R-R. In reality, 
migration contaminates a given radius R with stars coming from 
a variety of different initial radii R, however this measure can be 
interpreted as a way to quantify the intensity of radial migration 
for an initial (not null) gradient. If its value is greater than unity, 
for example (as it is the case in the inner disk in the phase of 
strong bar) then at any given radius, the metallicity variations in 
azimuth are similar to those expected from moving radially (in 
the center or anticenter direction) by more than 1 kpc. 



The ratio 6[Fe/H]/A[Fe/n] is also an indicator of the duration 
of the phase of strong migration: in the bar region, its value 
is indeed above unity for all the duration of high bar strength 
(cf Fig. [T] and [8]), thus indicating that inhomogeneities at any 
given distance from the center persist, i.e. stars are not mixed 
yet. The effect of mixing indeed starts at the end of the strong 
bar activity and migration phase (corresponding to t ~ 2 Gyr 
in our simulations): at this time, the ratio 6[Fe/H]/A\Fe/ti] dimin- 
ishes by a factor greater than 3, together with a decrease in bar 
and spiral arms strength. Once strong migration has ended, the 
timescale to disperse inhomogeneities of size AR in the disk is 
oc 2n/[(d£l/dR)AR], with D being the angular velocity of stars in 
the disk. For a fluctuation of size AR = 1 kpc, this corresponds 
to about 2 Gyr at R = 10 kpc, and 0.3 Gyr at R = 4 kpc. 

4. Conclusions 

In this work we study, by means of dissipationless, N-body sim- 
ulations, the signatures that radial migration induced by a bar 
imprints on the old stellar population of a disk galaxy. 

Firstly, we show that migration is initiated at the time of 
bar formation, and that it remains significant over the whole 
phase of high b ar str ength in agreement with previous results 
(iMinchev et afll201 lll20O:lBrunetti et alJl201 lb . 
Then we show that, as a consequence of this migration, signifi- 
cant azimuthal variations are generated in the metallicity distri- 
bution of old disk stars. 

At the peak of bar strength, the metallicity variations above or 
below the azimuthally averaged value are of the order of the ini- 
tial gradient. They are particularly strong in the bar region, with 
the higher metallicities observed along the bar major axis and 
the lower metallicities parallel to the bar minor axis. 
While the strength of the azimuthal variations depends on the 
initial stellar metallicity gradient, the ratio £[Fe/H]/A[F e /H] does 
not. This value is indeed independent on the initial metal distri- 
bution and, as we show, it is a measure of the strength of the bar 
and thus of the ongoing migration. Since migration, induced by 
bar and spiral arms, causes an inhomogeneous stellar distribu- 
tion in the stellar disk - the duration of this phase being quanti- 
fied by a ratio £[Fe/H]/A[F e /H] greater than unity - migration and 
mixing do not occur at the same time. Mixing, by definition, re- 
quires a homogeneous stellar distribution in the disk, and this is 
not compatible, as we show, with the phase of strong migration. 
Only when migration is over, mixing can be established. 
The presence of non-axisymmetric structures in a (face-on) 
galaxy disk does not provide information about the level of ac- 
tivity of the bar. It is not possible, for example, to determine, by 
means of morphological information alone, if a bar is in an early 
phase of evolution or if it has evolved already through a phase of 
vertical instabilities. Also, from morphology alone, is not possi- 
ble to quantify the level of stellar redistribution a disk is expe- 
riencing. Our work shows that metallicity variations in old stel- 
lar populatio ns of galaxy disks, possibly coupled with kinematic 
information (iGadotti & de Souzall2005l) . can be used as a diag- 
nostic tool for quantifying these processes, since their amplitude 
is strongly related to the level of activity of the bar and sub- 
sequent migration. Hence, chemical inhomogeneities over the 
whole disk would testify an ongoing strong bar activity, while 
inhomogeneities limited to the outer disk would suggest a fad- 
ing bar activity. Barred galaxies with axisymmetric distribution 
of metallicity could be the result of either a bar that has always 
been weak, or one that has now fade to weak activity, but note 
that the diff erence could still be detecta ble in the level of vertical 
dispersions (iGadotti & de Souz a 2005). 
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The predictions of our models can be detectable in barred 
galaxies with IFU surveys like CALIFA and ATLAS3D. In 
this sense, some encouraging results have been already pub- 
lished by means of long-slit s pectroscopy studies. Recently, 
ISanchez-Blazquez et al.l (120111) have indeed compared the 
metallicity properties of stars in the disk with those along the 
bar for two early-type galaxies, showing that in both cases bars 
show higher metallicities, and flatter gradients than disk stars in 
the same region. Our models predict the same trend and thus can 
be an important piece of information to add to the ongoing inter- 
pretation of these findings. 

Also in the Milky Way, azimuthal inhomogen eities have been re - 
ported, but currently only amon g young stars ([Luck et a"D(l2006l) . 
but see lLuck & Lambert! d2QlIl) : lDavies et al.l (l2009hT ~the varia- 
tions found may be thus the result of a patchy star formation. The 
detection (or not) of similar signatures in the old stellar compo- 
nent of the Galactic disk by future surveys like Gaia (and com- 
plementary spectroscopic data) may help in constraining the last 
epoch of strong activity of the bar and the amount of current ra- 
dial migration. 

Finally, since radial migration can be caused by several different 
physical mechanisms, as detailed in the introduction, it will be 
important to compare them to understand if all leave signatures 
similar to those discussed in this paper. 
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Appendix A: On the choice of initial conditions : are 
they realistic? 

Our initial conditions are compatible with a galaxy without 
prominent stellar asymmetries in its disk which experiences 
a phase of bar formation and strong bar activity, before the 
bar strength fades again. These models thus capture a possi- 
ble recurrent phase in the evolution of galaxy disks (see for ex 
Bournaud & Combes 2002). About the choice of adopting ini- 
tially a null azimuthal variation in the metallicity distribution, 
this is in agreement with the fact that at the end of the simula- 
tion, once the bar activity diminishes, azimuthal variations fade 
also, with typical values S ~ 0.4A, A being the metallicity gra- 
dient (see Fig and [8]). Our choice is thus compatible with our 
findings: the azimuthal variations fade, as the stellar distribution 
tends to axisymmetry. It can be noted also that the amplitude 
of azimuthal variations we find at the end of the simulation in 
the disk, for an initial A = -0.07dex/kpc is of the order of 2 
x 0.03 de*0= 0.06 dex (see Fig©. This value is lower than the 
dispersion induced in the system by blurring alone (see Fig. 1), 
thus ensuring that the choice of our initial conditions is realistic 
(any azimuthal variation of this order is included in the metal- 
licity dispersion at the pre-migration phase). To further prove 
that the choice of initializing the stellar disk with an axisym- 
metric metallicity distribution does not change our conclusions, 
we have post-processed our simulations, assigning, at any ra- 
dius, an initial azimuthal variation in metallicity S oc cos(20), 
i.e. z m (R,0) = ZolO~ omR lO s . The resulting metallicity variation 
maps are given in Fig. IA.ll for the case of 6 = 0.03cos(26), and 
show that even adopting initial azimuthal variations similar to 
those we get at the end of the simulation (S ~ 0.4A), we still re- 
cover the same trends found in Fig|7] The same results are found 
if the stellar disk has initially a more patchy metallicity distri- 
bution (Fig. IA.2b . Note that, for producing these maps, we have 
assigned non-zero azimuthal variations to the stellar disk at time 
t=0.4 Gyr. Indeed, any metallicity variation assigned at earlier 
times is completely washed out by shear before migration starts. 
Thus the amplitudes found in Fig. IA. 1 1 and IA.2I can be taken as 
upper limits. 



10 kpc 10 kpc 10 kpc 




Fig. A.l. Maps of azimuthal metallicity variations, £[F e /H] (dex), 
at different times for a stellar disk with azimuthal inhomo- 
geneities 6 ~ 0.4 A at t=0.4 Gyr (see text). Black contours cor- 
respond to equally spaced isosurfaces of the differential stellar 
density, See Fig|7]for comparison. 




Fig. A.2. Maps of azimuthal metallicity variations, S^q/u] (dex), 
at different times for a stellar disk with a patchy metallicity 
distribution at t=0.4 Gyr (see text). Black contours correspond 
to equally spaced isosurfaces of the differential stellar density, 
%diff- See Fig. IA. II and FigQfor comparison. 



6 The multiplication by a factor 2 is motivated by the fact that 5 cor- 
responds to half of the total variation 



